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We propose the mapping of polynomial of degree 25 constructed as a linear combination of 
powers of spin-S (for simplicity, we called as spin-S polynomial) onto spin-crossover state. The 
spin-S polynomial in general can be projected onto non-symmetric degenerated spin up (high-spin) 
and spin down (low-spin) momenta. The total number of mapping for each general spin-S 1 is given 
by 2(2 2S — 1). As an application of this mapping, we consider a general non-bilinear spin-S Ising 
model which can be transformed onto spin-crossover described by Wajnflasz model. Using a further 
transformation we obtain the partition function of the effective spin-1/2 Ising model, making a 
suitable mapping the non-symmetric contribution leads us to a spin-1/2 Ising model with a fixed 
external magnetic field, which in general cannot be solved exactly. However, for a particular case 
of non-bilinear spin-S Ising model could become equivalent to an exactly solvable Ising model. The 
transformed Ising model exhibits a residual entropy, then it should be understood also as a frustrated 
spin model, due to competing parameters coupling of the non-bilinear spin-S Ising model. 

PACS numbers: 05.50.+q; 05.70.-a; 64.60.My; 64.60.De 



I. INTRODUCTION 

One of the topics of great interest, in statistical 
physics and mathematical physics are the exact solv- 
able models. This is the case, i.e. for spin-1/2 Ising 
model without magnetic field solved at first in 1944 by 
OnsagerpJ, since that, the Ising model was widely in- 
vestigated using several approaches. On the other hand, 
higher order spin or even spin-1/2 Ising model with ex- 
ternal magnetic field are challenging issues in nowadays. 
Further exact solution were obtained only in a very lim- 
ited cases, mainly the honeycomb lattices [H, Some 
exact results has been obtained with restricted parame- 
ter, investigated by Mi and YangQ using a non-one-to- 
one transformation^. Therefore the non-bilinear spin- 
Slsing model that satisfy this transformation should ex- 
hibit frustrated states. 

The half-odd-integer spin Ising model already has 
been discussed previously by Tang[5j. Using the method 
proposed by Wujfl, Izmailian Q obtained an exact so- 
lution for a spin-3/2 lattice on square lattice with only 
nearest interaction or two body interaction spin, Izmail- 
ian and AnanikianQ also has been obtained an exact so- 
lution for a honeycomb lattice with spin-3/2. A particu- 
lar solution of these models could be obtained using the 
method proposed by JosephQ where any spin-5 1 could 
be projected onto a spin-1/2 Ising model. Another in- 
teresting method to map the spin-5 model onto spin-1/2 
Ising model has been proposed by HoriguchiJnJ. More 
recently we have obtained a set of rigorous mapping for 
half-odd-integer spin onto spin-l/2[ll|. where we have 
used a direct mapping, half of spin momenta are pro- 
jected onto spin down, while the remaining half spin 
momenta are projected onto spin up, we will call this 
process as spin-5 projection with symmetric degener- 
acy, this mapping is a non-one-to-one mapping. 

However, when we consider an integer spin we cannot 
perform the mapping symmetrically, this issue will be 
considered in this letter. Then using a non-symmetric 



projection we will discuss the mapping for a general 
spin-5 polynomial onto a spin-crossover state. 

On the other hand, the spin crossover (SC), some- 
times called as spin transition, is a phenomenon that 
occurs in some metal (i.e. Fe and Co) complexes wherein 
the spin state of the complex changes due to exter- 
nal perturbation such as a variation of temperature, 
pressure, light irradiation or an influence of a magnetic 
field[l3- The spin states of the atoms can change be- 
tween the high-spin (HS) state and low-spin (LS) state 
as a result of external stimuli, which can be understand 
as a non-symmetric degeneracy between HS and LS 
states. In this sense we find a equivalence between spin- 
5 polynomial and spin-crossover state. Another inter- 
esting equivalence could be also to that the metastable 
structure of a charge transfer phase transition, discussed 
by Miyashita et al.[ljj, where static metastability exist 
in a study of the charge transfer transition in the mate- 
rial (nC3H 7 )4N[Fe n Fe m (dto)3] (dto=C 2 2 S 2 ). 

The outline of this report is as follow: In sec. 2 we 
present the mapping of spin-5 polynomial onto spin- 
crossover state, in sec. 3 we apply to a a non-bilinear 
spin-5 Ising model and its relation to the metastable 
state, while in sect 4 we discuss the exactly solvable 
case. Finally in sec. 5 we present our conclusions. 

II. THE SPIN-S POLYNOMIAL 
TRANSFORMATION ONTO SPIN-CROSSOVER 
STATE 

In order to show the equivalence between spin-5 poly- 
nomial and spin-crossover state. Let us start considering 
as an example the projection of spin-3/2 polynomial, as 
follow, 

Cm 2 ( s ) = a , m + ai jm s + a 2: mS 2 + a 3iTO s 3 , (1) 

where a.i >m with i = 0, . . . , 3, are the coefficients to be 
determined using the projection of spin-3/2 onto spin- 
1/2, whereas by m we mean the number of solutions or 
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Table I. The projection spin-2 onto a non-symmetric spin-1/2 or spin-crossover state. By <?(+l) we mean the degeneracy of 
the spin up (HS), whereas g{— 1) means the degeneracy of spin down (LS). 



projectionsl $C m ^ 2(2 3 — 1). The polynomial spin-3/2 
of eq.([T]) can be rewritten alternatively by 



where 
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To find the coefficients of the polynomial spin-3/2 
eq. (|T|) , we use the following equation, 



a m = V _1 P. 
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each vector P m are defined as a column vector of matrix 



P = 




(6) 



the elements of matrix P, means the projection spin- 
3/2 polynomial onto spin up (+1) or high-spin (HS) and 
spin down (-1) or low-spin (LS). Each vector P m (col- 
umn of matrix P) are the permutations of at least one 
spin up and one spin down, here the columns of matrix 



P only represents the non-equivalent configurations of 

spin projections, for this case we have 5 "representative" 

configurations. We should recover easily the remaining 

configuration taken into account the exchange of spin 

(-) (-) 
s, <Jm (s) -H> <Jm (— s) and the global inversion of the 

(-) (-) 
polynomial ovn (s) O — <Jm (s). Therefore the spin-3/2 

polynomial mapping onto spin-1/2 are expressed by 
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The mappings given in ©-© already were consid- 
ered in reference 0, [Tlj . which corresponds to symmet- 
ric degeneracy mapping (the first 3 column of vector P). 
The remaining solutions corresponds to non-symmetric 
degeneracy, i.e. three spin momenta are projected onto - 
1 (LS), whereas the remaining spin moment is projected 
onto +1 (HS), or vice- verse. These solutions have not 
been considered yet in the literature. 

The next transformation that we discuss could be the 
spin-2 polynomial onto spin-1/2. It is not possible to 
map onto a spin-1/2 with symmetric degeneracy be- 
cause, we have five magnetic momenta to be mapped 
onto two eigenvalues ±1. Then the only possibility 
is to map by means of non-symmetric projection onto 
spin-1/2, this kind of mapping lead us to 30 poly- 
nomials. However we only need to obtain 9 "repre- 
sentative" polynomials, which are tabulated in table 
1. Once again using the exchange of magnetization 

(2) (2") 

<Jm (s) O a m (— s) and the global inversion of the poly- 
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nomial (Xm(s) O —am'(s), we could obtain easily the 
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In general the projection of spin-S polynomial onto 
cr(s) with non-symmetric spin degeneracy, we assume 
the following spin-S polynomial. 
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(s) 

where the coefficients Oj of the polynomial will be 
determined after projecting onto spin-1/2. 

To perform the spin-S polynomial projection we con- 
sider the Vandermonde matrix V' s ' with equidistant 
nodes [-S, S], whose elements of the node are Xj which 
corresponds just to the magnetic momenta of the spin- 
s', the elements of the matrix could be expressed ap- 



propriately as Xj 



-S + j, with j = 0,1, 2,..., 25, 



the explicit representation of the Vandermonde matrix 
is given by, 



The number of projection of spin-S polynomials that 
we can obtain are given by the permutations of the el- 
ements of the vector (|15[>. for each r = {1, ...,25}. 
Therefore, the total number of solutions for each spin-S 
is given by f^ 1 ) + f^ 1 ) + ■ ■ ■ + fff) = 2(2 2S - 1), 
thus m= {1,...,2(2 2S - 1)}. 

Using the matrix notation, we are able to write in 
general the spin-S polynomial, 



rL s ) = s&a& = s< s > f P^ 



(17) 



The inverse of the matrix \/( s > could be solved using 
the recursive equation presented recently by Eisinberg 
et aZ. |T3| , where was discussed a generic algorithm to 
obtain the elements of the inverse of Vandermonde ma- 
trix . Therefore the elements of the matrix \/^ are 



rewritten conveniently as in reference [13| . which reads 
as 
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we also define the vector aj? to represent the coeffi- 
cients for the all possible spin-S polynomial, while the 
elements of the vector P^ represents the projection 
with non-symmetric degeneracy HS and LS, which can 
be expressed by 
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where [ \ ] represent the first kind of Stirling num- 
ber, whereas F^^~ fe represents the hyper-geometric 
functionpj)]. 

Using the elements of inverse matrix v t 
to write the coefficient for each polynomial a(s), 



r(S) 



we are able 



3=0 



(19) 



note that here off^ are the elements of vector a.m\ 

The non-symmetric degeneracy projection of spin-S 
polynomial onto a spin-1/2 is given by 
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When the spin-S is half-odd-integer, we could recover 
the solution already obtained in reference [llj as a par- 
ticular case of our result, when the degeneracy becomes 
symmetric. 



by V we mean any permutation of the elements of P^', 
with r projections onto spin up (HS) and 25+ 1 — r pro- 
jections onto spin down (LS), assuming r = {1, . . . , 2S}. 
For a given r we have ( 2 r ) permutations, which projects 
the spin-S polynomial onto spin-1/2. We can verify that 
all matrix defined above have 2S + 1 dimension. 

In order to project the spin-S polynomial onto ±1 val- 
ues with non-symmetric degeneracy, we can use the ma- 
trix notation, so, the following algebraic system equa- 
tion becomes, 

PW=VW«W. (16) 



III. THE NON-BILINEAR SPIN-S ISING 
MODEL MAPPING ONTO WAJNFLASZ 
MODEL 

As an application of this mapping, let us consider 
the non-bilinear spin-S Ising model with two-body and 
high-order interaction term, whose Hamiltonian for ar- 
bitrary spin-S, can be written as 

25 2S 2S 

**= £ EEW'^-EE 5 ^ ( 21 ) 

<i,j> ki — 1 k2 — l i fc=l 
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with -ftTfci,fc 2 being the non-bilinear interaction terms, 
while by Bk corresponds to the high order anisotropy 
coupling. The < i,j > means the summation over the 
pairs of nearest-neighbor sites. 

In order to discuss the equivalence between non- 
bilinear spin-iS 1 Ising model and the Wajnflasz model(l8|. 
let us describe the Ising-like model or Wajnflasz 
modelH. 



III.l. The Ising-like model 

In order to describe the spin-crossover transition, we 
can use the Wajnflasz modeljl^, [2(| [2l|, where this 
model take into account the HS state and LS state, 
modeled simply by a nearest-neighbor interaction be- 
tween sites. 

Hence the Hamiltonian with arbitrary spin-S' Ising 
model can be mapped onto an effective spin-1/2 Ising- 
like model, which read as 



?C(W) 



E 

<i,j> 



Ja m (si)a m (sj) - h}^a m (si), (22) 



where J is the effective spin-crossover interaction pa- 
rameter, and whereas h corresponds to the effective ex- 
ternal magnetic field or the energy difference between 
HS and LS states. 

With aim of the eqs. (f2T|) and (f!2"2"]) becomes equiva- 
lent, we need to impose the following condition Us = 
'H l m {{<y}) — £q- Where the parameters must satisfy the 
relation below 



E' Q =MJal tm -Nhao !m , (23) 
B k =(h - jJa 0t m)a k ,mi k > 1, (24) 
Kk u k 2 =Jctk 1 , m ctk 2 ,m, ki > 1 and k 2 ^ I, (25) 

with M being the total number of nearest-neighbor spin 
pairs and TV being the total number of sites, while 7 
corresponds to the coordination number of the lattice. 

To study thermodynamics properties, we have to 
compute the partition function of Ising-like model, 



Z{fi) = exp(-/W s ), 



(26) 



where /3 = l/fegT, with ks being the Boltzmann con- 
stant and T the absolute temperature. 

Similar to that was discussed by Mi and YangQ, the 
eg. (f2"6"| can be rewritten as follow, 

Z m (J3) =e l3£ ° ^ g r (<7m,l)9r(<7m,2) ■ ■ ■ 9r(cTm,N)X 
{<x,}=±l 

eM-m' m {W})) (27) 
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,E I ln(^(^,,))-/3W^(W) 
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where H'miicr}) is the Hamiltonian of the effective spin- 
1/2 Ising model with non-symmetric degeneracy given 
by eq.(|22)l. 

For the purpose of mapping onto an exactly solvable 
model, we prefer to change the eq. ([2"5|) onto an usual 
standard form through a further transformation. There- 
fore, the partition function (f2"5|) becomes 

Z m (/3)=e^o J2 eM-P K^), (29) 
K}=±i 



<i,j> 



where 



~ ~p ( ln (9r(o- m ,i)) 4- ln(g r (a mJ ))) . (30) 

A further transformation, could leads us to an effec- 
tive Ising model Hamiltonian with temperature depen- 
dent field. 



III.2. The Ising model 

Finally using an additional transformation, the eq. 
([50)1 will be transformed onto a standard spin-1/2 Ising 
model, given simply by 



at this stage, r represents a standard spin-1/2, with J, 
h r and E r being parameters to be determined. Assum- 
ing the eqs. fBTJI) and (|31[) are equivalents, we have the 
following algebraic equations 
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—J - -5- (ln(ffr(l)) + MSr(-l))) = - J + E r . (34) 

P7 

Solving this algebraic system equations, we obtain the 
following relation 



J=J, 

h r = h - 



2(3 



In 
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Er = -— ln[.g r (l) 5r (-l)] 
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It is worth to highlight that when mapping is sym- 
metric, we have the following relation h r = h, similar to 
that discussed in reference fill ] - 

In addition the partition function of the Hamiltonian 
(f2"Tj) . is expressed by 



l-\/3; 
tanh(J*) = { y/2- 1; 

1/V3; 
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Some characteristic property of this model will 
be discussed now. At high temperatures, fc^T > 

2ft,/ In f g (i) ) i t ne term representing the effective field 



h r = h 



gr-(-l) 
9r(l) 



is positive and thus the spins 



have a positive expectation values (Tj) > 0. Whereas 
at low temperatures, fc^T < 2/i/ln^ 
negative expectation values (t, 



gr-(-l) 

9.(1) 



we have 



< 0. 



Miyashita et al.pjl discussed also the equivalence 
between spin-crossover phase transition and that the 
metastable structure of a charge transfer phase tran- 
sition, where static metastability exist in a study 
of the charge transfer transition in the material 
(nC 3 H 7 )4N[Fe n Fe m (dto)3] (dto=C 2 2 S 2 ). 



IV. THE EXACTLY SOLVABLE MODEL 



In order to that eq. ([38|) becomes an exactly solvable 
model, we consider the two-dimensional Ising modelQ, 
so, this model can be solved exactly when h r = 0, 
therefore it is equivalent to fix the magnetic field h = 

the 



27s ln 



gr-(-i) 

9.(1) 



in eq. (J22J). It is worth to note that, 

magnetic field h only depends of the degeneracy of spin 
momenta and is proportional to the temperature. 

The free energy of non-bilinear spin-S" Ising model in 
thermodynamic limit (N — > oo and M = jN/2), can 
be written in terms of standard spin- 1/2 Ising model by 
the following relation 



fm,S = -55 ln [ffr(-l) 



1— do,! 



9r{l) 



1+CKO,! 



+ /1/2, (40) 



where /jm means the free energy of spin- 1/2 Ising 
model. 

Therefore, for two-dimensional case, we can consider 
three types of lattice: triangular (IB), square [l| and 
honeycomb[l5| Ising model. 

The critical point (J*, h*) = (J/T c ,h r /T c ) for two- 
dimensional Ising model, in units of critical temperature 
T c , are given by 



and h* = 0, for triangular, square and honeycomb lat- 
tice, respectively. 

The critical points for honeycomb (7 = 3), square 
(7 = 4), and triangle (7 = 6) lattice can be fully re- 
covered, which is consistent with the results previously 
obtained by Mi and YangQ, for the case of spin-1 Ising 
model (for detail see table I of reference Q). However 
our result is quite general and is valid for any spin-5 
and for any coordination number. 

The non-bilinear spin-S* Ising model critical points 
should satisfy the relation below 



Bl 



(l ln (^w) ~ lJ*ao,m) Ofc,m, (42) 



(43) 



by * we mean the parameters are in units of critical 
temperature T c . 

On the other hand, the term E r and £q are responsible 
for the appearance of residual entropy, in other words 
this means due to competing parameters coupling of 
the non-bilinear spin-S* model could be considered as a 
frustrated spin model, therefore the entropy is given by 



|ln[5r(-l) 



1 — C*0,i 



1 + C*0,i 



(44) 



Using the above result (|4"4")l , we can obtain a residual 
entropy for the spin-1 Ising model discussed by Mi and 
Yang[4(, the first model has residual entropy given by 
S = ln(2) while the second model has no residual en- 
tropy. It is worth to notice the frustration properties of 
those model was not discussed by Mi and YangQ. 

Another simple example that we consider is the spin- 
3/2 Ising model, for the particular case such that satisfy 
the eas.([7HTT1). The residual entropy for eqs. ([7]) and 
([8]) are zero, inasmuch as the independent coefficients of 
those polynomials are ao,i = ao,2 = 0. However, for the 
last three polynomials we have a residual entropy given 
by S = ln(2), S = g ln(3) and S = ^ ln(3) for eqs. 
(|9l fTT|l respectively. It is interesting to highlight that, 
the residual entropy is independent of the coordination 
number or some other lattice structure parameters. 



V. CONCLUSION 

Different to those other methods developed to obtain 
this kind of results using a more involved approach, we 
have used a simple spin-S 1 polynomial projection onto 
spin-crossover state, with non-symmetric degeneracy of 
spin up or high-spin (HS) and spin down or low-spin 
(LS). In general the present projection obtained have 
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not be necessarily symmetric with relation to their spins 
up or down. Only as particular case of our results, 
we have the symmetric mapping which was previously 
considered in reference (Tl| . some additional results are 
found also using the decoration transformation method 
satisfying the 8- vertex model in our recent paper [l6j. 
Therefore we conclude that, there is a spin-S" polyno- 
mial transformation onto spin-crossover state, whose to- 
tal possible number of projection is given by 2(2 2S — 1). 
Through a further transformation we can map also onto 
a standard spin-1/2 Ising model. 

As an application of this mapping we consider the 
non-bilinear spin-S 1 Ising model which can be trans- 
formed onto spin-crossover state described by Wanjflasz 
model[l8j. Using a further transformation we obtain 
the partition function of the effective frustrated spin- 



1/2 Ising modelflj, making a suitable mapping this non- 
symmetric contribution leads us to a spin-1/2 model 
with a fixed external magnetic field temperature depen- 
dent given by eq. (|36[) . Therefore we conclude that, 
the non-bilinear spin-S 1 Ising model such that satisfy 
the projection proposed, must become equivalent to the 
Wanjflasz model[18|], with quite interesting properties 
such as residual entropy of the model, independent of 
the lattice structure. 
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